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摘 要 : 滩涂 是 我 国 重要 的 后 备耕 地 资源 ,为 探究 围垦 开发 对 沿海 滩涂 土壤 环境 的 影响 ， 对 浙江 省 慈溪 市 不 同 围 明 年 限 多 种 
土地 利用 方式 下 的 276 个 土壤 /沉积 物 样品 试验 检测 ,结合 Norm 值 优化 主 成 分 分 析 ， 筛选 出 有 机 质 (OM)、 和 氧化 钙 (CaO)、 硫 
(S). (Br), 9W(B). Pb), $&(Cd)7 项 指标 的 最 小 数据 集 ， 采 用 极 差 标准 化 和 变异 系数 定 权 的 Topsis 法 评价 土壤 质量 ,为 
揭示 阻碍 研究 区 土壤 质量 提升 的 障碍 因子 提供 参考 。 结 果 显 示 : 276 个 土壤 /沉积 物 样本 土壤 质量 综合 指数 为 44.6~74.1， 受 围 
星 时 间 影 响 差异 显著 , 呈 前 20a 逐 步 提高 一 20~30 a 相对 平稳 一 30~50 a 急剧 提高 的 趋势 ,从 利用 方式 来 看 ,土壤 质 量 呈 菜 地 > 
果园 > 农田 > 林地 > 养殖 > 潮 滩 > 荒地 的 趋势 。 农 业 活动 导致 的 增 肥 脱 钙 脱 盐 是 围垦 区 土质 提升 的 根本 原因 ， 但 研究 区 土壤 质量 
整体 仍 处 于 中 度 贫 壮 水 平 , 一 级 土壤 仅 占 研究 区 总 面积 的 7.1%， 质量 较 差 的 四 、 五 级 土壤 面积 超 一 半 。 较 低 的 OM 和 较 高 的 
CaO, Cd 是 限制 垦区 土质 提升 的 主要 障碍 因子 ,因此 ,土壤 培 肥 时 , 合理 施用 化 肥 农 药 是 进一步 提升 滩涂 围垦 区 土壤 质量 的 
关键 。 
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Soil quality evaluation in typical coastal reclamation zones based on weighted 
Topsis method" 


SUN Beiting, GAO Chao, ZHANG Yan" 
(College of Geographic and Oceanographic Sciences, Nanjing University, Nanjing 210046, China) 


Abstract: Tidal flat is an important reserve of cultivated land resource in China. There is high spatial heterogeneity, properties and 
environment changes of soil in the transformation from wetland to agricultural soil. In order to evaluate the effect of soil quality 
under different reclamation durations and multiple land use in a typical coastal area in Cixi County, Zhejiang Province, 276 
soil/sediment samples from farmland, vegetable land, orchard, forestland, aquafarm, wasteland and tidal flat reclaimed for 0 to 50 
years were collected and analyzed. Considering the spatial heterogeneity, principal component analysis was used in combination with 
Norm Value, then Minimum Data Set (MDS) with seven indexes was put forward. This included organic matter (OM), calcium oxide 
(CaO), sulfur (S), bromine (Br), boron (B), plumbum (Pb) and cadmium (Cd) under six reclamation duration and seven utility 
patterns. About 29.296 of the index accounted for 80.396 of the information. The Topsis method, which combined Maximum 
Difference Normalization method and Variable Coefficient Fixing Weight method, was used to assess the soil quality. Comparison 
with the Traditional Membership Function method, the used method a much more effective one. Besides, the distance to optimal 
vector directly explained the obstruction factor, which showed the primary obstacle was applicable in supervising soil environment 
and improving soil quality by the local government. The results suggested that the range of soil quality comprehensive indexes for 
the 276 soil/sediment samples was 44.6—74.1. The trend in the indexes increased sharply in first 20 years and then stabilized in the 
other 20-30 years, before it again increased for 30—50 years of reclamation. In general, the overall soil quality in the study area was 
moderate to infertile. The class I soil quality was only 7.196 of the samples. More than half of the soil samples belonged to classes IV 
or V. Based on land use types, soil quality was ranked in the following order: vegetable field > orchard > farmland > forest > 
aquafarm > tidal flat > wasteland. Paddy field was better than dry land for the same reclamation period, which indicated that basic 
paddy farmland was much more environmental friendly in Cixi reclamation district. Irrigation-induced decalcification and 
desalination effectively improved the quality of reclaimed soil. Meanwhile, low OM and high CaO and Cd primarily limited soil 
quality. Next were salt content as Br, S and B, followed by heavy metal as Pb. Soil content of Cd increased quickly to critical point of 
pollution, which was easily enriched in paddy and caused diseases. Reasonable use of chemical fertilizers was a key to further 
promotion of soil quality in reclaimed tidal areas. 
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滩涂 围 坚 区 是 受 海陆 交互 影响 的 半 人 工 半 自然 生态 系统 ， 是 我 国 重要 的 后 备 土地 资源 ,， 具有 很 强 的 开 
发 利用 潜力 。 围 垦 在 缓解 人 口 压力 、 保 持 耕 地 动态 平衡 等 方面 起 着 重要 作用 (4 E 50 a, 我 国 滩涂 围 尾 总 
面积 达 1.1x104~1.2x104km?, 成 为 滨海 围垦 力度 最 大 的 国家 中。 土壤 质量 是 “土壤 实现 各 种 功能 的 能 力 ”t6], 
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农学 上 侧重 维持 自然 植物 生长 的 能 力 I, 评价 土壤 质量 是 研究 农业 可 持续 发 展 的 重点 之 一 。 滨 海滩 涂 土壤 
具有 高 盐 低 养 、 地 下 水 矿 化 程度 高 等 特点 , 不 适宜 种 植 农作物 ,需要 在 围垦 前 期 实施 脱盐 补 肥 措 施 人 为 提 
高 土壤 质量 81。 这 种 人 为 干扰 极 大 地 改变 着 垦区 环境 。 晨 区 土壤 质量 偏 低 制 约 了 土地 的 开发 利用 中 ,滩涂 围 
晨 之 后 的 土壤 质量 演变 倍 受 关注 B101。 
受 土 壤 空 间 变 异性 影响 ,评价 土壤 质量 时 ,参评 因子 的 取舍 至 关 重 要 。Andrews 等 UI 和 Yemefack 55U?1 
首次 提出 以 主 成 分 分 析 中 最 高 因子 荷载 值 选取 最 小 数据 集 (MDS) 作 为 评价 最 小 单元 , 广 受 学 者 认可 。 李 桂 
MUZE MDS 基础 上 引入 矢量 常 模 (Norm 值 )， 既 能 最 大 程度 代表 所 选 参数 又 尽 可 能 减少 信息 损失 , 但 
指标 仅 包 括 土壤 物理 性 质 和 养分 , 涵盖 信息 较 单 一 。 贡 赖 等 ( 习 、 邓 绍 欢 等 9 引入 生物 指标 丰富 了 最 小 数据 
集 , 但 未 全 面 涉及 土壤 基本 理化 性 质 。 岳 西 杰 等 7 依托 GIS 县 加 多 种 地 形 部 位 、 植 被 覆盖 和 施肥 等 构建 了 
土壤 质量 评价 MDS， 却 忽视 了 时 间 因 素 的 影响 。 已 有 研究 建立 的 最 小 数据 集 多 集中 于 土壤 元 素 和 地 形 地 貌 
因素 , 而 参评 因子 随 围垦 时 间 和 土地 利用 的 变化 与 轨迹 等 方面 的 研究 其 少 。 

评价 土壤 质量 多 基于 隶属 函数 B10, 但 隶属 函数 中 指标 同一 化 转折 点 的 取 值 均 由 专家 打分 确定 ,存在 
明显 的 主观 性 ; 且 在 计算 土壤 质量 综合 指数 时 ， 多 采用 层次 分 析 法 或 主 成 分 法 来 确定 权重 ， 其 中 两 两 判断 
和 矩阵 的 建立 由 专家 确定 ， 主 成 分 定 权 通常 只 取 蓓 载 较 高 的 若干 因素 ， 因此, 定 权 也 具有 一 定 主观 性 。 已 有 研 
究 对 于 内 部 参数 的 选择 仅 通 过 公式 换算 稍 加 完善 , 鲜 有 学 者 引入 新 的 模型 进行 评价 。 

本 文选 用 综合 荷载 Norm 值 并 引入 浙江 省 慈溪 市 海 涂 围垦 区 土地 利用 和 年 限 变化 计算 综合 得 分 ， 筛 选 
出 评价 最 小 单元 MDS, 结合 滩涂 围 愤 区 受 自然 演 苦 和 人 为 活动 影响 强烈 的 特点 ,丰富 指标 体系 ,以 达到 婚 
最 大 程度 代表 所 选 参数 又 尽 可 能 减少 信息 损失 的 目的 ， 并 采用 变异 系数 加 权 Topsis 方法 ， 从 时 空 二 维 角 度 
客观 评价 围垦 区 土壤 质量 ,揭示 阻碍 研究 区 土壤 质量 提升 的 障碍 因子 。 
1 材料 与 方法 
11 研究 区 概况 

浙江 省 慈溪 市 共 修筑 11 条 围垦 海 是 0819。 本 文 研 究 明 区 位 于 慈溪 市 东南 部 (32"07~32?19'N， 
121922'121?34'E), HAR 114 km?, 包含 5 条 主要 围垦 海 堤 ,该 区 属 北 亚热带 南 缘 季风 季候 , 温和 湿润 ， 光 
热 充足 , 年 均 气 温 16.3 C, 年 均 日 照 1 933.5 h, 年 均 降 水 量 1 325 mm。 主 要 土壤 类 型 为 滨海 盐 土 、 灰 潮 土 
和 水 稻 士 ， 母 质 为 滨海 相 沉积 物 。 研 究 区 地 理 位 置 及 土地 利用 情况 见 图 1。 
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图 1 研究 区 地 理 位 置 及 土地 利用 情况 


Fig. 1 Geographical location and land use types of the study area 


12 数据 获取 

时 外 采样 于 2013 年 7 月 完成 , 按 研究 区 面积 通过 网 格 法 均匀 布点 ， 共 采集 276 个 土壤 /沉积 物 样 品 。 采 
样 时 使 用 手持 GPS 定位 , 在 目标 点 10 m 范围 内 , 五 点 法 采集 0~20 cm 的 表层 土壤 /沉积 物 , 混合 后 形成 一 
个 土壤 样品 。 土 壤 样 品 经 自然 风干 后 过 10 目 得 储存 备用 。 根 据 遥 感 影像 上 不 同时 期 的 海 堤 海岸 线 确 定 围 
RIK, 按 围 明 年 限 自 海洋 向 陆地 方向 分 别 标记 为 IT 区 、I 区 、JII 区 、IV 区 、V 区 、VI 区 ; 按照 土地 利用 
方式 , 样品 取 自 农田 111 个 、 菜 地 47 个 、 果 园 48 个 、 林 地 17 个 、 养 殖 27 个 、 匾 地 21 个 、 淹 滩 5 个 (图 
NE 1)。 


R1 研究 区 土地 利用 及 采样 点 分 布 情况 
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Table 1 Land use and distribution of soil samples in the study area 


ERAS PAMO ” 海 坦 修 筑 年 从 HEER 
District Numbers of Year for building Years of reclamation 
number samples dams (a) 

I 21 RELÆ Tidal flat 0 
II 49 2002—2012 10 
TIT 25 1990—2002 20 
IV 13 1980—1990 30 
V 35 1970—1980 40 
VI 133 —1970 50 


土壤 样品 理化 测试 项 目 包括 粒 径 比 ( 粉 砂 粒 / 茜 粒 )、 pH、 有 机 质 (OM)、N、P、K20、CaO、Na;0、MegO、 
Al>03、Fex03、SiO0，、Cl、S、Mn、Br、B、Cu、Zn、Pb、Cd、As、Hg 和 Cr。 采 用 多 目标 区 域 地 球 化 学 
调查 规定 的 分 析 方 法 及 质量 监控 方案 20， 由 江苏 省 地 质 调查 研究 院 测 试 应 用 研究 所 、 安 徽 省 地 质 实验 研究 
所 承担 分 析 测 试 。 
13 土壤 质量 评价 体系 建立 

根据 慈溪 市 滩涂 围 导 区 土壤 理化 性 质 (养分 、 常量 与 微量 元 素 、 土壤 质地 等 ), 通过 主 成 分 分 析 (PCA 法 ) 


T 提取 因子 荷载 ,计算 矢量 常 模 (Norm 值 ) 通过 多 变量 方差 分 析 (MANOVA)、 均 值 多 重 比较 (S-N-K 比较 法 ) 
三 ” 等 方法 结合 围 明年 限 和 士 地 利用 方式 建立 围垦 区 土壤 质量 的 最 小 数据 集 。 

e 综合 荷载 Norm 值 可 弥补 PCA 中 单一 化 导致 的 信息 缺失 。Norm 值 的 几何 意义 为 变量 在 多 维 空间 中 矢 
N 量 常 模 的 长 度 , Norm 值 越 大 , 表明 变量 对 所 有 主 成 分 的 综合 荷载 越 大 01, 涵盖 信息 就 越 综合 。 

Lo k 

e Na = JS xA) (1) 

N 1 

KC 0 式 中 : Ni 为 第 i 个 变量 的 Norm 值 , 指 在 特征 值 1>1 TITEL E ERA er wR i MERES k 
v 0 个 主 成 分 上 的 荷载 ; 4 为 第 大 个 主 成 分 的 特征 值 。 

eo 采用 多 变量 方差 分 析 (MANOVA) 衡 量 土 地 利用 变化 对 各 土壤 参数 的 影响 ,调整 后 的 决定 系数 R? 表征 
CN 利用 方式 对 各 参数 的 影响 大 小 ,均值 多 重 比 较 (S-N-K 比较 法 ) 估 算 围垦 年 限 对 各 土壤 参数 的 影响 用 最 早产 
二 生 显 著 影响 的 年 限 了 作为 时 间 对 研究 区 土壤 理化 性 质 产 生 影响 的 界限 09。 

S< 为 使 最 终 因 子 得 分 具有 可 比 性 , 对 Norm 值 ,R2、 了 进行 简单 线性 变换 ， 其 中 Norm 值 和 R? 为 正 指标 , 其 


© 值 越 大 表明 对 土壤 参数 的 影响 越 大 ， 因 此 ,最 大 值 赋 为 1， 其 他 变量 除 以 相应 的 最 大 值 ; 7 为 逆 指 标 , 进行 相 
= 反 转 化 。 求 和 得 到 每 个 土壤 参 数 的 总 分 值 。 
特征 值 4>1 的 各 个 主 成 分 中 , 若菜 因子 荷载 >0.5 则 归 为 该 主 成 分 组 ,若菜 参数 同时 在 两 个 主 成 分 中 载 
© 荷 >0.5, 或 在 每 一 个 主 成 分 中 荷载 都 <0.3， 则 单独 归 为 新 的 一 组 。 为 进一步 减少 数据 元 余 , 抽取 分 组 后 在 最 
高 得 分 10% 范 围 内 的 参数 ,观察 其 与 小 组 间 其 他 参数 的 相关 性 (x)， 若 高 度 相关 (xr>0.5) 且 在 0.01 置信 水 平 下 
E, 则 只 选取 得 分 最 高 的 进入 MDS, 反之 则 全 部 进入 MDS. 
1.4 土壤 质量 评价 方法 
Topsis 法 是 一 种 逼近 理想 解 的 多 属性 排序 法 20,， 它 根据 比较 样本 与 最 佳 解 和 最 劣 解 的 相对 距离 来 评价 
样本 的 优 劣 。 该 方法 对 样本 容量 无 特殊 要 求 ， 且 不 受 参考 序列 选择 的 干扰 , 已 广泛 应 用 于 水 质 综 合 评价 和 
农业 综合 开发 项 目 评价 R1351, 但 在 土壤 质量 评价 中 应 用 甚 少 。Topsis 法 主要 步骤 为 : 
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COE 2 b VES 
20 C min(x;] 
正 指标 : ” max[ xX; }—min{ x; } (2) 
max{ xXj}- X; 
à, 三 i 
逆 指 标 : ^o max(x,)—-min(x;] (3) 
| max[ x;] - x; -x| 
ij 
适中 指标 : max[x,] m 


式 中 : max{xjy} 、min{xy}、xo; 分 别 为 j 指 标 中 的 最 大 值 、 最 小 值 、 适 中 值 。 
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根据 标准 化 矩阵 A 所 有 样本 实际 数据 ,由 各 指标 的 最 优 值 与 最 劣 值 确定 最 优 样本 4+ 与 最 劣 样本 4-。 
最 优 样本 : A+ =(ai,at,...,ai) (5) 
最 劣 样本 : A =(oia ar) (0) 
式 中 : a+、4a 分 别 为 正 、 负 理想 解 。 
加 变异 系数 定 权 


(7) 
式 中 : 为 第 j 个 指标 的 变异 系数 , 5) 为 第 j ARREA, ox 7958 j 个 指标 的 平均 值 。 
@ 计 算 各 土壤 样本 与 最 优 及 最 劣 样本 的 加 权 欧 氏 距 离 D* 与 D- 


Dr = yo (da (8) 


式 中 : D* 5j D- 分 别 表示 第 i 个 评价 对 象 与 最 优 方案 及 最 劣 方 案 的 距离 , aug 表示 第 i 个 评价 对 象 在 第 j 个 评价 
指标 上 经 极 差 标准 化 后 的 取 值 ，%; 表 示 第 j 个 指标 的 权重 。 
计算 土壤 质量 综合 指数 Ci 


D; 


C; kx 
D} + D; 


(9) 
式 中 : Ci; 表示 第 i 个 样本 土壤 质量 综合 指数 ， 即 每 个 土壤 样本 与 最 优 土壤 样本 的 接近 程度 ;k W100. C; BUR 
便 在 0 到 100 之 间 , C; 越 接近 100, 说 明 该 土壤 样本 越 接 近 最 优 土壤 质量 , 土壤 质量 越 高 ， 反之, 土壤 质量 越 
差 。 
1.5 ”障碍 因子 诊断 模型 

评价 土壤 质量 可 直观 反映 围 殿 对 土壤 的 改 展 效果， 有 针对 性 地 消除 阻碍 土质 提高 的 障碍 因子 对 促进 围 


坚 区 土壤 发 育 有 重要 意义 。 本 文 用 Topsis 评价 公式 (8) 中 “与 最 优 方案 的 距离 D** 即 可 表示 第 i 个 土壤 样本 的 
障碍 度 ， 几 何 意义 直观 , D! 越 大 则 障碍 度 越 高 ， 越 需 对 该 样本 进行 改良 。 
2 结果 与 分 析 
2.1 围垦 区 土壤 基本 属性 

研究 区 土壤 总 体 呈 弱 碱 性 , pH 均值 8.3, 随 围 虹 时 间 呈 先 上 升 之 后 陆续 下 降 的 趋势 ,变异 系数 (CV)3.8%， 
表明 在 人 工 围 娠 后 土壤 pH 变化 不 大 。Al203、Fe203、SiO; 的 平均 含量 分 别 为 128.90 gkg, 50.39 gkg 
636.89 gkg, 符合 弱 碱 性 土壤 特性 。 选取 粉 砂 粒 与 颖 粒 的 比值 作为 土壤 机 械 组 成 参数 031， 其 变异 系数 高 达 
55.1196, 差异 较 大 ,表明 研究 区 土壤 机 械 组 成 受 人 为 活动 影响 显著 。 参 考 土壤 比重 换算 方法 29， 土壤 容重 
均值 为 1.474 gcm-3,， 略 高 于 全 国 大 多 数 农用 地 区 , 但 仍 属于 结构 较 好 的 土壤 , 适宜 发 展 农业 。 有 机 质 平 均 
含量 13.4 gkg!, 参照 全 国 第 2 次 土壤 普查 土壤 肥力 分 级 275， 肥 力 状况 属 中 低 等 水 平 。 

Cl 受 海水 影响 极 大 ， 有 具有 高 度 变异 性 (CV>200%)， 且 含量 与 围垦 时 间 高 度 相 关 ， 治 海 堤 逐 步 变 化 ,说 
明 以 农业 为 主 的 人 类 活动 加 速 了 滩涂 沉积 物 脱盐 过 程 。 受 海 相 动力 作用 主导 , S 和 Br 含量 高 于 全 国平 均 水 
平 ， 由 海 向 陆 存 在 较 大 的 差异 。 

WAX Cu, Zn, Pd, Cd, As, Hg, Cr 七 大 主要 重金 属 含量 略 高 于 全 国土 壤 背 景 值 ， 虽 未 达到 污染 水 
平 , (HH Cu. Pd. Cd. As. Hg 变异 系数 超 15%。 
2.20 最 小 数据 集 选 取 

构建 合理 的 指标 体系 是 Topsis 法 的 关键 在 主 成 分 分 析 (PCA) 基 础 上 ， 爱 加 多 元 方差 分 析 (MANOVA)、 
均值 多 重 比较 (S-N-K 比较 法 ) 和 Spearman 相关 分 析 等 ， 得 到 最 小 数据 集 。 

主 成 分 特征 值 4>1 的 4 项 主 成 分 累积 贡献 率 高 达 80.29%, 可 解释 绝 大 部 分 指标 的 变异 性 ， 按 每 个 参数 
在 各 主 成 分 的 载荷 分 成 5 组 。 运 用 公式 (1) 计 算 各 参数 的 Norm 值 , 经 线性 变化 后 求 和 得 到 综合 分 值 ( 表 2)。 

表 2 各 参数 综合 得 分 及 土壤 质量 评价 最 小 数据 集 确定 
Table 2 Scoring of soil parameters and determining of Minimum Data Set for soil quality assessment 
Xa m 初始 分 值 mitial value 线性 变换 Linear transformation 综合 得 分 
Norm 值 调整 R2 Norm 值 调整 RR 


Soil parameter Group T T Scoring 
Norm value Adjusted R? Norm value Adjusted R? 


m 


a 
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Taste 
Siluclay 
KO 
AlO3 
Fe;0s 


1 2.136 


1 3.092 
1 3.120 
1 3.120 
1 3.020 
1 2.800 
1 2.805 
1 2.314 
1 2.700 
3.022 
1.869 
1.992 
2.040 
2.289 
1.758 
1.731 
1.022 
2.319 
2.326 
2.481 
2.702 
2.942 
1.967 


Un aa ta ta ta ta DD NNNN 一 


2.123 


0.158 4 


0.291 
0.338 
0.332 
0.292 
0.189 
0.183 
0.057 
0.283 
0.332 
0.516 
0.458 
0.743 
0.477 
0.223 
0.457 
0.1 


N UU QN — tn d d tA t — QN ta tA tA wm tA 


0.333 
0.338 
0.388 
0.585 
0.466 
0.721 


Un N A tn U 上 上 一 


0.543 


0.685 


0.991 
1.000 
1.000 
0.968 
0.897 
0.899 
0.741 
0.865 
0.968 
0.599 
0.638 
0.654 
0.734 
0.563 
0.555 
0.327 
0.743 
0.745 
0.795 
0.866 
0.943 
0.630 
0.680 


0.213 


0.391 
0.455 
0.447 
0.393 
0.254 
0.246 
0.077 
0.381 
0.447 
0.695 
0.616 
1.000 
0.642 
0.300 
0.615 
0.135 
0.448 
0.455 
0.522 
0.787 
0.627 
0.970 
0.731 


0.25 


0.20 
0.33 
0.20 
0.20 
0.20 
0.17 
1.00 
0.20 
0.20 
0.25 
0.25 
0.20 
1.00 
0.17 
0.33 
0.50 
1.00 
0.25 
0.33 
0.20 
0.25 
0.50 
0.20 


1.147 


1.583 
1.788 
1.647 
1.561 
1.352 
1.312 
1.818 
1.446 
1.615 
1.543 
1.505 
1.854 
2.376 
1.030 
1.503 
0.962 
2.191 
1.450 
1.651 
1.853 
1.820 
2.101 
1.611 


1 组 中 Pb 综合 得 分 最 高 ， 且 与 组 内 其 他 参数 相关 度 高 ， 因 子 得 分 前 10% 的 


AbOs(r-0.832"). Fe20;3(r=0.815*)5 Pb 高 度 相 关 , 结合 Norm 值 大 小 需 剔 除 元 余 ， 故 第 1 组 中 只 选取 Pb 


进入 MDS; 第 2 组 类 似 ， 只 选取 Cd 2 


组 中 综合 得 分 最 高 ， 


低 ， 故 在 选取 OM 的 同时 , 添加 MgO 和 Br 进入 MDS, 


碱 过 程 中 Ca 与 Mg 之 间 有 很 高 的 相似 性 09, 
最 终 进 入 研究 区 土壤 质量 评价 最 小 数据 集 的 有 效 参数 为 Pb、Cd、S、 


经 围垦 年 限 差异 性 显著 检验 ,MDS 7 项 参数 的 组 内 变异 均 小 于 纪 
质 的 影响 远大 于 内 部 随机 误差 。 结 合 MDS 中 不 同 


体 呈 先 下 降 后 上 升 


BULL E EHE E ER 


含量 整体 呈 先 上 升 后 下 降 又 上 升 再 下 降 的 波动 变化 趋势， 
43.4 gkg! 后 速 降 至 围垦 50 a 的 29.4 g.kg- ， 脱 钙 脱 碱 高 月 


35 


30 


N 
[] 


有 机 质 OM (g kg-1) 
S 


. 

. s~ 

中 

. . e 
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. [| 

15 [] xd 
5 i e 
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与 前 10968] Cl(r=-0.515””) 高 度 相关 ， 


但 与 MgO(/=-0.162") 和 Br(/=0.343”*) 的 相关 度 较 
H MgO 存在 较 弱 的 变异 性 (CV<10%)， 盐 渍 土壤 脱 


因此 用 Cao &HX MgO SE TIG ES] ES DX ERES B AD. i E PS, 


proT] 
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HEER 


30 
ERR 
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Reclamation years 


Tr - - 
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围垦 年 限 
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ER 
Reclamation years 


2) HERAN 


B. OM. Cao 和 Br 等 7 项 指标 。 
有 间 变 异 , 说 明 人 为 围 展 对 土壤 理化 性 
围垦 年 限 土壤 化 学 性 质 特征 (图 


含量 整 


r4 H 


的 变化 趋势 ， 滩涂 沉积 物 有 机 质 平均 含量 为 12:0 g-.kg !, BRER, HERE, LRA 
E 10 a 后 的 8.8 gkg, 之 后 逐步 增 至 围垦 50 a 后 的 15.2 gkg! HIX RZ. Cao 平均 
由 沉积 物 中 的 40.8 gkg 382] EFEM] E& 30 a 的 
现象 显著 ,适宜 作物 生长 。 


bal ga 
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铅 Pb (mg kg-1) 
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图 2 


围垦 年 限 
Reclamation years 


最 小 数据 集 参 数 在 不 同 围垦 年 限 下 的 土壤 化 学 性 质 特 征 图 


Fig.2 Chemical characteristics of soil under different reclamation years in the Minimum Data Set 


常 微量 元 素 S 和 Br 的 变化 相似 , 在 围 虹 前 后 呈 先 下 后 升 再 下 的 “S 型 ?波动 , 组 间 差 异 明显 , 均值 分 别 
从 滩涂 时 期 的 0.44 gkg! Fi 18.52 mg kg"! KEH] 50 a 的 0.21 g-kg-! 41 7.22 mg-kg-!, 减少 了 近 一 半 。B 


的 变化 恰 相 反 , Sis S 型 "波动 , 但 总 体 平稳 , E 


E it Ji M 66.81 mg: kg"! 小幅 升 至 68.76 mg.kg-1。 


5| E& 10 a 后 ， 重金属 元 素 Pb 平均 含量 从 滩涂 沉积 物 时 32.6 mg.kg- 下 降 至 22.2 mg.kg 1， 减少 31.9%， 


之 后 的 40 a 间 略 升 至 26.7 mg:kg-!。Cd 虽 也 在 围垦 10 a 出 现 最 低 值 0.12 mg- kg, 但 在 之 后 40 a 间 波 动 上 
升 至 0.20 mg.kg-1， 比 潮 滩 时 期 增加 36.2%, 临近 土壤 环境 质量 (GB15618 一 2008) 一 级 标准 28]。 


23 土壤 质量 评价 


A MDS 的 7 项 参数 按照 土壤 理化 性 质 可 分 成 肥力 指标 (OM)、 常 量 及 微量 元 素 (CaO0、S、 了 Br、 


B) 和 重金 属 指标 Pb、Cd)( 表 3). Hb, OM 是 正 指 标 ， 其 值 越 大 ,土壤 质量 越 高 ; Pb 与 Cd 为 逆 指 标 , 其 值 


越 高 , 污染 越 严重 ,土壤 质量 越 差 ; 其余 指 标 均 为 适度 指标 , 适中 值 的 确定 在 研究 区 元 素 理化 性 质 的 基础 


可 


上 综合 考虑 包括 中 长 江平 原 区 (GBW07429-GSS15)P9、 土 壤 成 分 分 析 元 素 含量 参考 值 (GSS1-16)P2、 浙 江 省 


和 全 国土 壤 元 素 背景 值 80 在 内 的 国家 标准 。 
通过 Microsoft Excel 2013 3 


ÍT Topsis 评价 ， 先 对 MDS 中 的 7 项 指标 按 式 (C2)-(4) 进 行 极 差 标准 化 ,转化 
至 0-1 范围 内 , 使 评价 模型 低 敏感 高 稳健 便于 操作 。 
K, 选用 变异 系数 按 式 (8) 定 权 既 能 保持 客观 公正 ， 


围 晨 活动 受 时 间 影响 强烈 ,不 同 海 堤 内 元 素 含 量变 化 
又 能 增强 分 区 内 指标 的 可 比 性 。 


表 3 研究 区 土壤 质量 评价 体系 


Table 3 Evaluation system of soil quality in the study area 


标 层 准则 层 指标 层 权重 属性 适中 值 
Target layer Criterion layer Factor layer Weight Attribute Modest value (mg:kg"!) 
(o GEHURR 肥力 Fertlty OM —— 038 —— E Posie 0/1 

Soil quality 常 、 微 量 CaO 0.194 适中 Modest 1530 

Macro-element and S 0.082 适中 Modest 176 

mcro-element Br 0.097 适中 Modest 6.25 
B 0.059 适中 Modest 63 
重金 属 Pb 0.070 XÉ Inverse / 
Heavy mental Cd 0.160 XÉ Inverse / 
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采用 式 (5)-(9) 计 算 每 个 样本 土壤 质量 综合 评价 值 Ci, 值 越 高 表示 土壤 质量 越 高 。 研 究 区 276 个 样本 点 


的 土壤 质量 综合 评价 值 在 44.6 至 74.1 之 间 , 平均 值 56.1, 最 优 土壤 质量 样本 是 种 植 水 稳 (Oryza sativa) 的 农 


H, 而 最 劣 土壤 质量 


样本 为 跳 跳 鱼 养殖 塘 。 


土壤 质量 整体 变化 趋势 受 年 限 主导 , 局 部 变化 特征 受 利用 方式 影响 (图 3)。 研究 区 土壤 质量 与 年 限 高 度 
相关 (R=0.98), 50 a 间 提 高 20%， 呈 “前 20 a 逐步 提高 一 20~30 a 相对 平稳 一 30~50 a 急剧 提高 ”的 轨迹 变化 。 
从 利用 方式 来 看 ,土壤 质量 菜 地 > 果园 > 农田 > 林地 > 养殖 > 潮 滩 > 范 地 。 农 田 、 菜 地 和 果园 的 土壤 质量 综合 评 


价 均值 达 57.6， 由 于 人 为 施肥 量 大 ,利用 效益 明显 高 于 养殖 、 荡 地 和 潮 滩 ,说 明 开发 农业 对 滩涂 土壤 脱 碱 增 


肥 有 较 好 的 改良 作用 。 
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图 3 不 同 围垦 年 限 和 利用 方式 下 研究 区 土壤 质量 综合 评价 值 


Fig.3 Comprehensive indexes of soil quality under different reclamation years and land-use types in the study area 


步 按 作物 细 


分 为 种 植 瓜 果 蔬菜 [西瓜 (Citrullus lanatus), $8?) (Vitis vinifera)、 桃 子 (Amygdalus persica)、 


南瓜 (Cucurbita moschata)、 毛 豆 (Glycine max)、 油 菜 (Brassica campestris). 3ESE(Allium ordorum) 等 ] 的 旱地 


和 种 植 水 稻 的 水 田 两 种 。 以 既 有 水 田 又 有 旱地 的 开 区 为 例 , 在 时 间 固 定 条 件 下 , 水 稻田 平均 士 质 (C=52.1) 
优 于 旱地 平均 士 质 (C=51.9),， 恰 与 研究 区 整体 水 田 劣 于 旱地 的 结果 相反 (图 3)， 且 在 区 区 出 现 类 似 情况 ， 表 
明 在 同一 时 间 内 种 植 水 稻 的 环境 效益 更 优 。 因 此 ,当地 农民 在 种 植 瓜 果 以 追求 经 济 利益 的 情况 下 ,应 重视 
保护 环境 , 合理 施用 化 肥 农 药 , 避免 人 为 输入 过 量 的 重金 属 元 素 Cd 诱发 地 方 性 骨 痛 病 等 。 

林地 土壤 质量 呈现 较 高 均值 是 因为 种 植 的 万 年 青 (Tlex chinensis) (Melia azedarach) 等 绿化 树 相 比 


mug 


THIBDNETESE GA, 不 需要 过 多 的 人 为 照料 , 肥力 不 及 农用 地 ,重金 属 含量 也 相对 少 。 土壤 质量 潮 滩 略 高 于 


mihi, 说 明 围垦 后 的 


土地 应 充分 合理 利用 ,避免 空地 总 废 , 政府 应 出 台 鼓 励 及 补助 政策 减少 土地 资源 的 滔 


费 。 


~ 


在 全 国 范 围 内 比 


N 
人 土地 利用 Land use 
a 
* 
. 
* W Tidal flat 
土壤 质量 综合 评价 值 Ci 
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4 研究 区 土壤 质量 分 布 情况 
Fig.4 Distribution of soil quality in the study area 


较 , 研究 区 土壤 质量 仍 呈 中 度 贫 将 状态 1。 对 样本 点 插值 后 参考 ArcGIS 中 自然 断 点 
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法 (enks),， 将 土壤 质量 等 级 分 成 5 种 (图 4: 一 级 土壤 最 少 , 仅 8.09 km?， 占 总 面积 的 7.1%， 大 多 分 布 在 围垦 
时 间 最 长 的 V 区 ,以 农田 和 菜 地 为 主 ; 二 级 土壤 占 19.596, 广泛 分 布 于 最 早 的 3 REREN, MUR R 
bd. 林地 为 主 , 辅 之 以 少量 农田 和 养殖 用 地 ; 三 级 土壤 占 总 面积 的 20.7%, 主要 是 农田 ， 散 布 于 II 区 和 VI 区 ; 
占 27.4% 的 四 级 土壤 和 25.2% 的 五 级 土壤 分 布 在 最 新 修筑 的 两 条 海 担 内 ， 为 少量 农田 、 林 地 加 上 绝 大 多 数 
养殖 、 荒 地 和 淹 浴 。 农 用 地 的 土壤 质量 远 高 于 养殖 和 荒地 ， 说 明 近 50 a 间 的 围 海造 田 不 仅 缓解 了 粮食 和 人 
口 的 压力 ,带动 了 当地 的 经 济 ， 且 对 土壤 环境 以 及 促进 生态 系统 的 发 育 都 有 积极 影响 。 

2.4 土壤 质量 障碍 因子 诊断 

基于 Topsis 改进 后 的 障碍 因子 诊断 模型 ， 计 算 MDS 中 土壤 各 项 指标 的 障碍 度 。 研 究 区 整体 因子 障碍 
程度 排序 为 OM>CaO>Cd>Br>S>B>Pb, 结合 表 4 的 分 区 情况 可 知 ,在 不 同年 限 和 利用 方式 下 ,主要 障碍 因 
子 不 尽 相 同 。 总 体 来 看 OM 对 土壤 质量 的 阻碍 占 主导 ; 海 相 沉 积 物 受 OM. Cao 和 Br 影响 较 大 ,围垦 前 
10 a 呈现 出 短暂 的 低 肥 高 钙 变 化 ,OM 和 Cao 阻碍 加 剧 , 土壤 逐渐 脱离 海水 环境 , Br 含量 下 降 不 再 成 为 主要 
障碍 ; 围垦 20-40 a, Cd 的 负面 作用 持续 加 剧 , 成 为 阻碍 土质 的 主要 因素 ; 50 a 殿 区 土壤 完成 脱盐 脱 钙 达 到 
稳定 状态 ,肥力 持续 提高 , Cd 成 为 仅 次 于 OM 的 第 2 障碍 因子 ， 比 上 一 个 10a 提高 32.9%， 是 未 围 必 时 期 的 
2.9 倍 , 需 引 起 重视 。 苔 加 土地 利用 方式 来 看 , 旱地 虽 提 高 了 肥力 但 Cd 含量 也 随 之 增加 。 


表 4 研究 区 土壤 质量 指标 障碍 度 


Table 4 Limiting degree of soil quality indexes in the study area 


土地 利用 方式 Land use type OM CaO S Br B Pb Cd 
滩涂 Tidal flat 0.163 0.056 0.01 0.022 0.001 0.003 0.009 
围垦 10 a Reclamation for 10 a 0.223 0.061 0.002 0.003 0.002 0.001 0.005 
围垦 20 a Reclamation for 20 a 0.166 0.057 0.002 0.002 0.002 0.001 0.011 
FE] E& 30 a Reclamation for 30 a 0.139 0.068 0.008 0.001 0.001 0.001 0.010 
围垦 40 a Reclamation for 40 a 0.124 0.045 0.001 0.004 0.001 0.001 0.017 
围垦 50 a Reclamation for 50 a 0.115 0.019 0.001 0.002 0.002 0.002 0.026 
水 田 Paddy field 0.188 0.06 0.004 0.005 0.002 0.001 0.008 
旱地 Dry land 0.117 0.027 0.001 0.003 0.001 0.002 0.025 
整体 Total 0.145 0.038 0.002 0.004 0.002 0.001 0.018 


3 讨论 与 结论 

采用 主 成 分 分 析 、 多 元 方差 分 析 等 方法 , 结合 Norm 值 、 围 明年 限 和 利用 方式 ， 从 24 项 土壤 常规 理化 
性 质 中 筛选 出 OM. CaO, S. Br. B. Pb. Cd 7 项 指标 作为 衡量 研究 区 土壤 质量 的 最 小 数据 集 MDS。7 
项 指标 全 面 履 盖 肥 力 、 常 微量 元 素 和 重金 属 三 大 方面 , 仅 用 29.2% 的 指标 就 包含 了 80.3% 的 信息 ， 既 减少 了 
工作 量 又 降低 了 测试 成 本 ,为 后 续 评价 围 暴 区 土壤 质量 提供 了 一 个 简洁 全 面 的 模型 框架 。 

滩涂 围垦 活动 下 ,土壤 性 质 受 人 为 影响 发 生 明 显 变化 。OM 含量 先 下 降 后 上 升 的 变化 趋势 , CaO 含量 在 
30 a 急剧 下 降 后 至 50 a 趋 于 稳定 ， 这 与 已 有 研究 相符 。 晨 区 整体 仍 属 轻 度 贫 凌 状态 , 说 明 围垦 后 开发 需 经 
历 较 长 的 脱 碱 增 肥 阶段 ， 实现 滩涂 向 农田 的 转换 还 需 一 定 的 时 间 积 累 和 人 为 改良 , Cd 含量 在 开发 后 期 快速 
上 升 甚至 接近 污染 临界 点 ,， 这 是 已 有 探究 围 明 区 土壤 所 忽略 的 一 点 ， 易 富 集 于 水 稻 的 Cd 元 素 需 引起 广泛 
的 重视 ,应 从 种 植 区 的 选择 到 撤 药 施肥 等 阶段 全 面 监测 农作物 中 Cd 含量 。 

相 比 于 传统 基于 隶属 函数 计算 加 权 关 联 度 的 诊断 模型 ， 变异 系数 定 权 的 Topsis 法 减少 主观 打分 且 几 何 
意义 更 加 直观 。 研 究 区 土壤 质量 与 围垦 年 限 高 度 相关 ， 呈 “前 20 a 逐步 提高 一 20 一 30 a 相对 平稳 一 30 一 50 a 
急剧 提高 ”的 趋势 .从 土地 利用 方式 来 看 , 土壤 质量 为 菜 地 > 果园 > 农田 > 林地 > 养殖 > 淹 滩 > 荒地 ， 说 明 开 发 农 
业 对 滩涂 土壤 脱 钙 增 肥 有 明显 改良 作用 。 同 一 围垦 年 限 下 水 田 土质 优 于 旱地 ,表明 荔 溪 滩涂 更 适宜 种 植 水 
稳 以 保护 基本 农田 建设 。 农 业 活 动 导致 的 增 肥 脱 钙 脱 盐 是 围 惰 区 土质 提升 的 根本 原因 , 但 研究 区 整体 土质 
仍 属 于 中 度 贫 竣 水 平 , 一 级 土壤 最 少 , 仅 占 总 面积 的 7.1%, 质量 最 差 的 四 、 五 级 土壤 超 一 半 。MDS 中 , K 
子 障碍 程度 排序 为 OM>CaO>Cd>Br>S>B>Pb, 较 低 OM 含量 和 较 高 的 CaO 、Cd 含量 是 暴 区 土质 提升 的 主 
要 限制 因子 。 淡 水 浇灌 是 初期 土壤 质量 提升 的 根本 原因 ， 养 分 提高 以 及 脱 钙 加 剧 是 后 期 继续 提升 的 关键 ， 
fH 50 a 间 Cd 含量 增加 了 近 3 倍 。 因 此 在 增加 土壤 培 肥 的 同时 , 合理 施用 化 肥 农 药 是 进一步 提升 滩涂 围垦 
区 土壤 质量 的 有 效 措施 。 


土壤 最 小 数据 集 因子 的 选 定 目 前 仍 处 于 摸索 阶段 ， 本 文 指标 选取 时 未 包含 生物 指标 , 也 未 考虑 农业 生 
产 过 程 (如 灌溉 与 施肥 ) 的 影响 。 随 年 限 增加 和 土地 利用 变化 土壤 的 空间 变异 性 将 加 强 , 因此， 建立 一 个 全 
面 评价 围 晨 区 土壤 质量 演变 的 指标 体系 ,以 及 采样 密度 对 模型 稳健 性 的 影响 还 有 待 后 期 深入 研究 。 
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